Exciton spin-flip rate in quantum dots determined by a modified local density of 

optical states 
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The spin-flip rate that couples dark and bright excitons in self-assembled quantum dots is ob- 
tained from time-resolved spontaneous emission measurements in a modified local density of optical 
states. Employing this technique, we can separate effects due to non-radiative recombination and 
unambiguously record the spin-flip rate. The dependence of the spin-flip rate on emission energy 
is compared in detail to a recent model from the literature, where the spin flip is due to the com- 
bined action of short-range exchange interaction and acoustic phonons. We furthermore observe a 
surprising enhancement of the spin-flip rate close to a semiconductor-air interface, which illustrates 
the important role of interfaces for quantum dot based nanophotonic structures. Our work is an 
important step towards a full understanding of the complex dynamics of quantum dots in nanopho- 
tonic structures, such as photonic crystals, and dark excitons are potentially useful for long-lived 
coherent storage applications. 

PACS numbers: 78.67.Hc, 78.47.-l-p, 72.25.Fe 



Understanding the optical properties of soHd-state 
quantum emitters is much required due to their increas- 
ing importance for all-solid-state quantum photonic de- 
vices for quantum information processing [Ij]. Spin de- 
grees of freedom of electrons and holes in semiconduc- 
tor quantum dots (QDs) impose exciton fine structure: 
The long-lived dark-exciton states can recombine after a 
phonon-mediated spin-flip process, whereby a radiating 
bright exciton is formed [2, l3| . Dark exciton recombina- 
tion is essential to explore: it influences the performance 
of QD single-photon sources 0, Q and the quantum- 
optical properties of QDs [6*1 , and must be accounted for 
when assessing inhibition of spontaneous emission in pho- 
tonic crystals 0, [1] ■ Furthermore, long-lived spin states 
in the solid state are required for spin qubits and 
dark excitons might be attractive in this context since 
they can be prepared and read out optically. 

In a QD, the exciton spin lifetime is greatly extended 
compared to a bulk medium or a quantum well due to 
the discrete energy spectrum. The mechanism behind 
the spin-flip processes has been debated in the literature 
[loL , and experimental tests for unbiased QDs have 
been lacking. Here we employ controlled modifications 
of the local density of optical states (LDOS) to deter- 
mine the dark exciton spin-flip rate from time-resolved 
spontaneous emission measurements. The LDOS is mod- 
ified by placing QDs at controlled distances from a 
semiconductor-air interface, which was previously used 
to distinguish radiative and non-radiative recombination 
processes of bright excitons Controlling the LDOS 
provides a powerful tool to obtain detailed insight into 
QD dynamics; thus cavity quantum electrodynamics was 
recently employed for QD spectroscopy in the regime of 
continuous pumping [l3| . We observe a surprising depen- 
dence of the spin-flip rate on the distance to the sample 



surface ranging several hundred nanometers, and a sensi- 
tive dependence on the QD emission energy. The latter is 
compared in detail to theory describing the dependence 
of the spin-flip rate on QD size, and provides valuable 
insight on the exciton spin-flip mechanism in QDs. 

Figure [Ija) illustrates the fine structure of the low- 
est exciton state for InAs/GaAs QDs The exci- 
ton is formed from the conduction band electron state 
(spin 1/2) and the heavy-hole valence band state (to- 
tal angular momentum 3/2). As a result, four exciton 
states are formed. They are characterized by the pro- 
jections of the total angular momentum onto the growth 
axis, which attain the values ±1 or ±2 for bright and 
dark excitons, respectively. The splitting between dark 
and bright energy levels Adb is determined by the ex- 
change coupling between electron and hole spins, and 
the bright-exciton levels are typically a few hundred |j,eV 
above the dark-exciton levels. The bright excitons \b) de- 
cay to the ground state \g) (no excitons) by emission of 
a photon through a dipole-allowed transition with a rate 
7rad- Radiative transitions from dark excitons \d) are for- 
bidden, however, they can decay through a phonon me- 
diated spin-fiip process (rate 7db) transforming the dark 
exciton into a bright exciton. The rate of the reverse pro- 
cess is denoted 7bd- Finally non-radiative recombination 
was recently proven to be significant for self-assembled 
QDs I4I, i.e., bright and dark excitons can recombine 
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re- 



non-radiatively with rates denoted 
spectively. 

We have measured spontaneous emission decay dynam- 
ics of InAs QDs positioned at 28 different distances from 
a GaAs/air interface, see inset of Fig. [T](b). Further ex- 
perimental details can be found in Refs. [S]- We 
record the number of photons emitted per time interval 
N{t) from the QDs. Solving the rate equations for the 
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Figure 1; (Color online) (a) Three-level exciton scheme con- 
sisting of a bright \b) and a dark |d) state that are coupled 
through the spin-flip rates 7bd and 7db- Radiative recombi- 
nation (7rad) to the ground state \g) (no exciton) is only pos- 
sible for the bright state, while both the bright (7nrad) and 
the dark state (7nrad) can decay non-radiatively. (b) Typical 
decay curve acquired at a detection energy of 1.204 eV (blue 
squares) at the sample with QDs positioned 170 nm from the 
interface. The red line is a fit of the bi-exponential model to 
the data. Inset Sketch of the sample showing three different 
distances between the QDs and the GaAs-air interface, (c) 
The weighted residual obtained from the bi-exponential fit to 
the data resulting in Xr ~ 1-11. 



three-level system of Fig. [Ija) results in a bi-exponential 
decay N{t) = Afe"'^'* + Ase~'^''*. Here subscripts s and 
f denote slow and fast decay rates, respectively. Fig- 
ure [ijb) shows an example of a spontaneous emission 
decay curve together with a bi-exponential fit obtained 
with a fixed background level to account for the measured 
contribution from dark-counts and after-pulsing of the 
detector. Excellent agreement between the experiment 
and the model is apparent from the weighted residuals 
reproduced in Fig.[T]^c). 

The four parameters 7f/s,^f/s are obtained by fitting 
the experimental decay curves. We have 7f ~ 7rad + 7nrad 
and 7s ~ 7nrad' where it has been assumed that the spin- 
flip rate is slow compared to both radiative and non- 
radiative processes, which is a very good approximation 
as seen below. Thus, the spin-flip rate cannot be obtained 
from the measured rates 7f /s since they are dominated by 
other decay contributions. However, the amplitudes Affs 
contain additional information and the dark-bright spin- 
flip rate is contained in their ratio 



Aijz) ^ 7f(z)-7sPbft^0) 
As{z)~ jdhiz) Pd{t = 0) 



(1) 



where pb(^ — 0)/Pd{t — 0) is the ratio between the 
initial populations of bright and dark excitons after an 
excitation pulse and z is the distance to the interface. 
The repetitive nature of the experiment implies that 
the slow amplitude contains contributions from previous 
excitation pulses, which is accounted for by correcting 
the amplitudes according to As As (l ~ e^^"'^) where 
T = 13.5 ns is the repetition period of the excitation 
laser. The spin-flip rate is obtained from experimentally 



ldh{z) = 



7f(^)-7s Pb(t = 0) 
1 + Af{z)/As{z) Pd(t = 0)- 



(2) 



It is essential that the ratio of the amplitudes enters in 
Eq. ^ since the absolute size of Af/g will depend on, e.g., 
the total collection efficiency of the radiated light or the 
number of QDs probed. Since the bright and dark ex- 
citons recombine radiatively on the same transition (cf. 
Fig. [TJa)) these dependencies do not contribute to the 
amplitude ratio. Quite remarkably, assessing the am- 
plitudes allows for extraction of spin-flip rates that are 
slower than the repetition period of the measurement. 

In Fig. [2] the slow decay rate and the ratio of the am- 
plitudes are plotted versus distance to the interface. The 
fast rate (not shown) varies periodically with distance in 
full agreement with the modified LDOS caused by the re- 
flecting interface [13] ■ In contrast, the slow decay rate is 
constant within the error-bars of the measurement, which 
confirms that it is dominated by non-radiative recombi- 
nation of dark excitons (7ii].ad)- average decay rate 
of 7s = 0.097 ± 0.008 ns^^ matches very well with the 
non-radiative decay rate of bright excitons at the same 
emission energy, which is expected since the bright and 
dark exciton binding energies are very close. Af{z)/As{z) 
is expected to vary proportional to the LDOS through 
7f(2;) (cf. Eq. ([T])) and indeed clear oscillations with a 
period matching the LDOS is observed. In addition the 
spin-flip rate is found to vary, which gives rise to the 
overall increase of the amplitude ratio with z. Assum- 
ing an exponential decrease of the spin-flip rate with z, 
which will be discussed further below, and multiplying 
by the calculated LDOS we can model the experimental 
data very well, see Fig. [D^b). 

In order to extract the spin-flip rate from Eq. we 
need to estimate the ratio of the initial populations of 
bright and dark excitons. Due to the non-resonant exci- 
tation, the feeding of the QDs is non-geminate [l3| , thus 
dark and bright excitons are formed with equal probabil- 
ity. A slight imbalance between the initial populations is 
nonetheless created due to the finite probability of creat- 
ing biexcitons. The recombination of biexcitons is dom- 
inated by radiative decay, which always leaves behind a 
bright exciton and thus increase pb(0) relative to pd(0). 
The long lifetime of the dark excitons moreover increases 
the probability of biexciton formation since a residual 
population of dark excitons is persistent when the sub- 
sequent excitation re-excites the QDs. In the present 
experiment the excitation density was fixed such that on 
average 0.1 excitons were generated per QD and by solv- 
ing the rate equations for the QD population including 
the biexciton level we find ph(t = 0)/pd{t = 0) ~ 1.25. 

The spin-flip rate versus distance to the interface is pre- 
sented in Fig. [3] for six different emission energies. It is 
found to vary between 4.7|js~"'^ and 14.7|j.s~^ correspond- 
ing to long-lived exciton spin lifetimes between 68 ns and 



3 



0.15 



^ 0.10 
(0 



O 0.05 
CO 

0.00 
150 



< 100 

50 



(a) 



(b) 



4* 



100 150 200 250 300 

Distance z [nm] 




1.18 1.20 1.22 1.24 1.26 1.28 

Emission energy [eV] 



i» »; «» * 



. xrm 
1* ♦ * «♦ ♦ 



50 



■ 1.272 eV • 
1.204 eV 


1.252 eV 
1.187 eV 


A 1.216 eV 
• 1.170 eV 




100 150 


200 


250 



300 



Distance z [nm] 



Figure 2: (Color online) (a) The slow decay rate as a function 
of distance to the interface measured at 1.204eV. No system- 
atic dependence on distance is observed confirming that the 
slow rate is due to non-radiative decay of dark excitons. The 
blue line indicates the averaged rate, (b) Measured ratio of 
the fast and slow amplitudes versus distance. The dashed 
blue curve is obtained by multiplying the calculated LDOS 
with an overall exponential decrease of the spin-flip rate with 
the distance to the interface. 



215 ns. Surprisingly the spin-flip rate is observed to de- 
pend exponentially on distance to the interface, which is 
observed for all emission energies. The characteristic de- 
cay length varies between 24nm and 106nm. This increase 
of the spin-flip rate in the vicinity of the interface could 
be caused by enhancement of acoustic phonons at the in- 
terface. In this case, we estimate the phonon wavelength 
from Aph = hv / Adh where h is Planck's constant, v is 
the sound speed and Ajb the exchange energy splitting 
between bright and dark states that must be matched by 
the phonons. For GaAs we estimate Aph 65 — 85 nm for 
longitudinal phonons using the range of exchange ener- 
gies relevant in the experiment, which matches the length 
scale observed in the experiment. 

The above results illustrate clearly the importance of 
nearby surfaces when seeking a quantitative understand- 
ing of the dynamics of QDs in nanostructures. Recent 
examples of intricate surface effects include the increased 
emission rate observed near a semiconductor interface 
[l5| . while charge trapping near surfaces was suggested 
as the mechanism responsible for the surprisingly large 
QD-nanocavity coupling efficiency at very large detun- 
ings observed under non- resonant excitation 17 , [isl . [l9| . 
These are examples of the new physics found in solid- 
state implementations of quantum electrodynamics ex- 
periments. 

Let us now consider the energy dependence of the spin- 



Figure 3: (Color online) The spin-flip rates versus distance 
z to the interface for six different emission energies. The 
plotted curves are flts to the experimental data assuming an 
exponential decay of the spin-flip rate with the distance to 
the interface. The inset shows the exponential decay length 
versus emission energy. 
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Figure 4: (Color online) Measured spin-flip rate at different 
emission energies for the sample with {z — 302 nm) (solid 
black circles). The theoretically predicted spin-ffip rate using 
the model of Ref. [13] for parameters that reproduce the 
energy dependence of the radiative decay rate (dashed blue 
line) and for optimized parameters (solid red line). 



flip rate, which allows to assess how the spin-flip process 
varies with QD size. In order to minimize the effect of 
the interface, we plot in Fig. 2] the spin-flip rate versus 
energy for the sample with largest QD-interface distance 
{z = 302 nm). The spin- flip rate is observed to increase 
significantly with the emission energy, and varies from 
6ns"i at 1.170 eV to 13 ns^^ at 1.272 eV meaning that 
the characteristic spin-flip time can be prolonged by more 
than a factor of 2 by varying the QD size. 
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The energy-dependence of the spin-flip rate can be 
compared to theory using the model introduced in Ref. 
[lo| . In this work the dark-bright exciton spin- flip rate 
was calculated due to the combined effect of short-range 
exchange interaction and acoustic phonons. A complex 
behavior is expected since the spin-flip rate is predicted 
to depend on the energy splitting between the lowest 
and first excited exciton states A12, the exchange en- 
ergy splitting Adb, as well as the size of the electron and 
hole wavefunctions. 

The experimental data for the energy-dependent spin- 
flip rates are compared to theory using experimentally 
realistic parameters: By recording the emission spectrum 
of the QDs at high pump power where the ground state 
is saturated, we obtain A12 = f 20 ± 20 meV. The ex- 
change splitting depends strongly on the indium mole 
fraction ^] and it is enhanced in QDs due to the strong 
localization of the electron and hole wavefunctions as de- 
scribed by the enhancement factor j2fj. Furthermore, 
the spin-flip rate depends on the electron and hole wave- 
functions through a form factor Q . We have calculated 
the electron and hole wavefunctions using the theory and 
parameters of Ref. [T3| . All other parameters are taken 
from Ref. [13] and the references given above. 

The comparison to theory is presented in Fig. 21 The 
dashed blue curve is for parameters that reproduce the 
frequency dependence of the radiative decay rate [3l (as- 
pect ratio: f/2, indium mole fraction: 0.46, QD heights: 
between 4 nm and 5.8 nm). The observed increase of 
the spin-flip rate with energy is reproduced by the the- 
ory. However, the calculated rates are two orders of 
magnitude smaller than the measured values. We have 
explored the origin of this deviation by optimizing the 
model parameters in order to maximize the spin-flip rate, 
see the solid red curve in Fig. S] (aspect ratio: f/3, in- 
dium mole fraction: 0.90, QD heights: between 1.9 nm 
and 2.4 nm). In this case the predicted spin-flip rates 
approach the measured rates, although a systematic dis- 
crepancy is still observed. Furthermore, a rather weak 
frequency dependence of the spin-flip rate is predicted, 
since the finite confinement potentials impose lower lim- 
its to the achievable compression of the wavefunctions. 
We conclude that additional spin-flip processes must be 
included in a complete theory. Additional contributions 
may arise due to spin-orbit coupling 11| or long-range 
exchange effects [23 1. 

In conclusion, we have measured the rate for spin flip 
from dark to bright excitons using time-resolved fluores- 
cence spectroscopy. The spin-flip rate increases signifi- 
cantly when approaching the sample surface suggesting 
an enhancement due to surface acoustic phonons. The 
energy dependence of the spin-flip rate was compared 
to a recent theory, where the spin flip is induced by 
short-range exchange interaction between electrons and 
holes while the required energy is provided by acoustic 
phonons. Our results illustrate the importance of tak- 



ing the QD fine structure into account when interpreting 
luminescence experiments, and will be important in or- 
der to obtain quantitative understanding of light-matter 
interaction in QD based optical devices. 

We gratefully acknowledge the Danish Research 
Agency for financial support (projects FNU 272-05-0083, 
272-06-0138 and FTP 274-07-0459). 



[1 
[2: 

[s: 

[4: 

[s: 

[6: 
[7: 

[8 

[9: 

[lo; 
[11 

[12 

[13 
[14 

[15 

[16 
[17 

[18 

[19 



* Electronic address: 'p elo@fotonik.dtu.dk. | 
B. Lounis and M. Orrit, Rep. Prog. Phys. 68, 1129 
(2005). 

O. Labeau, P. Tamarat, and B. Lounis, Phys. Rev. Lett. 
90, 257404 (2003). 

J. M. Smith, P. A. Dalgarno, R. J. Warburton, A. O. 
Govorov, K. Karrai, B. D. Gerardot, and P. M. Petroff, 
Phys. Rev. Lett. 94, 197402 (2005). 

A. Strauf, N. G. Stoltz, M. T. Rakher, L. A. Coldren, 
P. M. Petroff, and D. Bouwmeester, Nat. Photon. 1, 704 
(2007). 

T. Lund-Hansen, S. Stobbe, B. Julsgaard, H. 
Thyrrestrup, T. Siinner, M. Kamp, A. Forchel, and P. 
Lodahl, Phys. Rev. Lett. 101, 113903 (2008). 
M. Reischle, G. J. Beirne, R. Rofibach, M. Jetter, and P. 
Michler, Phys. Rev. Lett. 101, 146402 (2008). 
P. Lodahl, A.F. van Driel, I.S. Nikolaev, A. Irman, K. 
Overgaag, D. Vanmaekelbergh, and W.L. Vos, Nature 
430, 654 (2004). 

B. Julsgaard, J. Johansen, S. Stobbe, T. Stolberg-Rohr, 
T. Siinner, M. Kamp, A. Forchel, and P. Lodahl, Appl. 
Phys. Lett. 93, 094102 (2008). 

A. Imamoglu, D. D. Awschalom, G. Burkard, D. P. Di- 
Vincenzo, D. Loss, M. Sherwin, and A. Small, Phys. Rev. 
Lett. 83, 4204 (1999). 

K. Roszak, V. M. Axt, T. Kuhn, and P. Machnikowski, 
Phys. Rev. B 76, 195324 (2007); erratum: Phys. Rev. B 
77, 249905(E) (2008). 

E. Tsitsishvili, R. v. Baltz, and H. Kalt, Phys. Rev. B 
72, 155333 (2005). 

J. Johansen, S. Stobbe, I. S. Nikolaev, T. Lund-Hansen, 
P. T. Kristensen, J. M. Hvam, W. L. Vos, and P. Lodahl, 
Phys. Rev. B 77, 073303 (2008). 

M. Winger, A. Badolato, K.J. Hennessy, E.L. Hu, and A. 
Imamoglu, Phys. Rev. Lett. 101, 226808 (2008). 
M. Bayer, G. Ortner, O. Stern, A. Kuther, A. A. Gor- 
bunov, A. Forchel, P. Hawrylak, S. Fafard, K. Hinzer, 
T. L. Reinecke, S. N. Walck, J. P. Reithmaier, F. Klopf, 
and F. Schafer, Phys. Rev. B 65, 195315 (2002). 
S. Stobbe, J. Johansen, P.T. Kristensen, J.M. Hvam, 
and P. Lodahl, submitted, cond-mat ar Xiv:0902. 0344 1 /1 
(2009). 

B. Baylac, X. Marie, T. Amand, M. Brousseau, J. Bar- 
rau, and Y. Shekun, Surface Science 326, 161 (1995). 
K. Hennessy, A. Badolato, M. Winger, D. Gerace, M. 
Atatiire, S. Guide, S. Fait, E. L. Hu, and A. Imamoglu, 
Nature (London) 445, 896 (2007). 

M. Kaniber, A. Laucht, A. Neumann, J. M. Villas-Boas, 
M. Bichler, M.-C. Amann, and J. J. Finley, Phys. Rev. 
B 77, 161303(R) (2008). 

A. Laucht, F. Hofbauer, N. Hauke, J. Angele, S. Stobbe, 



5 



M. Kaniber, G. Bohm, P. Lodahl, M.-C. Amann, and J.J. 
Finley, New Journ. of Physics 11, 023034 (2009). 
[20] H. Fu, L.-W. Wang and A. Zunger, Phys. Rev. B. 59, 
5568 (1999). 

[21] R. Romestain and G. Fishman, Phys. Rev. B 49, 1774 
(1994). 

[22] A. Grodecka, L. Jacak, P. Machnikowski, and K. 



Roszak, in "Quantum Dots: Research Developments", 
edited by P. A. Ling (Nova Science, New York, 2005); 
|arXiv:cond-mat / 0404364 
[23] J.W. Luo, A. Franceschetti and A. Zunger, Phys. Rev. B 
79, 201301(R) (2009). 



